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CROSS-CORRELATION (C 2 ) IMAGING FOR 
WAVEGUIDE CHARACTERIZATION 
ROMAN A. BARANKOV 

ABSTRACT 

Confined geometries, such as optical waveguides, support a discrete set of eigen- 
modes. In multimoded structures, depending on the boundary conditions, superposi- 
tion states can propagate. Characterization of these states is a fundamental problem 
important in waveguide design and testing, especially for optical applications. 

In this work, I have developed a novel interferometric method that provides com- 
plete characterization of optical waveguide modes and their superposition states. The 
basic idea of the method is to study the interference of the beam radiated from an 
optical waveguide with an external reference beam, and detect different waveguide 
modes in the time-domain by changing the relative optical paths of the two beams. 

In particular, this method, called cross-correlation or C 2 -imaging, provides the 
relative amplitudes of the modes and their group delays. For every mode, one can 
determine the dispersion, intensity and phase distributions, and also local polarization 
properties. 

As a part of this work, I have developed the mathematical formalism of C 2 -imaging 
and built an experimental setup implementing the idea. I have carried out an exten- 
sive program of experiments, confirming the ability of the method to completely 
characterize waveguide properties. 
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Chapter 1 
Introduction 



In optical applications, the demand for reliable waveguide characterization methods 
comes from the development of several fiber-based optical technologies. First, the 
network traffic in optical fiber communication systems demands further increase in 
optical channels to overcome the limitations of currently deployed wave-division mul- 
tiplexing (WDM) systems, such as the bandwidth of optical amplifiers and input 
power. 



Mode-division multiplexing (MDM) pioneered in Ref. (Berdague and Facq, 1982) 



and recently developed in Refs.(Hanzawa et al., 2011; Ryf et al., 2011; Salsi et al., 



2011) is one possible solution of the network-capacity problem. MDM systems em- 



ploying several fiber modes for transmission of information require real-time monitor- 
ing of modal power decomposition. 

Several approaches have been suggested to realize the monitoring, such as the far- 



field imaging of the output facet of multimoded fibers (Rittich, 1985), and a method 



using ultrafast sources and specially designed probe fibers to temporally resolve the 



launched power (Golowich et al., 2004). A promising alternative for modal power 



monitoring in real-time is provided by the tilted-fiber Bragg gratings. The gratings 
couple guided modes to radiation modes encoding, among other characteristics, the 



modal power distribution in the probed systems (Wagener et al., 1997; Westbrook 



et al, 2000; Feder et al, 2003; Yang et al., 2005; Yan et al, 2012). 



High-performance, high-power fiber-lasers (Richardson et al., 2010) is another 



area in which waveguide characterization is of outmost importance. The design of 
novel fiber-platforms capable of achieving excellent output beam-qualities requires 



the propagation of only one mode in fibers that are not strictly single-moded (Dong 



et al., 2009| |Ramachandran et al., 2006} |Galvanauskas et al., 2008t |Stutzki et al., 



2011; Koplow et al., 2000). A usual measure of the laser-beam quality in these 
systems is the so-called M 2 -parameter (Siegman, 1990a; Siegman, 1990b). However, 
being essentially an integral characteristic, it gives a very rough estimate of higher- 



order mode content (Wielandy, 2007). 



A recently reported characterization technique - S imaging (Nicholson et al., 



2008; Nicholson et al., 2009) - enables the direct measurement of modal content 



of multimoded waveguides by recording spatially resolved output of the fiber in 
frequency domain. Historically, this approach appeared from the early studies of 



mode dispersion of higher-order modes (Menashe et al., 2001) and multi-path in- 



terference QRamachandran et al., 2003; Ramachandran et al., 2005; Ramachandran, 



2005). In S -imaging method, one of the modes, usually the fundamental one, is used 



as a "reference" mode for the analysis of the interference signal, which is a funda- 
mental impediment of the approach. In particular, in this case, a reasonable estimate 
of the relative power levels of the modes may be obtained only when the "reference" 
mode has the largest power in the propagating beam. Moreover, the presence of 
multiple interferences of every mode with all other modes complicates the analysis of 
the recorded signal and makes the power reconstruction unreliable. As a result, the 
method fails to characterize the generic case of multiple modes having similar power 
levels. Yet another difficulty is related to the inherent insensitivity of the method to 
the polarization properties of the modes, which introduces uncontrollable approxima- 
tions in the interpretation of the results, in the general case of elliptically polarized 
states. 
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An alternative approach of optical low-coherence interferometry does not suffer 
from the limitations of S 2 -imaging, since it employs the external reference beam with 
well-defined characteristics, which interferes with all the modes of the test waveg- 



uide (Nandi et al., 2009; Ma et al., 2009). This technique does not require any 
prior knowledge of group-delay or dispersive characteristics of the reference mode. 
In contrast to S 2 imaging, all the modes having arbitrary relative power levels and 
polarization properties can be measured independently of one another. 

In this work, I have developed and realized this idea. In particular, I have worked 
out the mathematical formalism forming the basis of the method and demonstrated 
that interferometric signal contains complete information about waveguide modes, 
including modal decomposition, intensity and phase distributions of all the modes, 
and also their polarization properties. The formalism was applied for characterization 
of several optical waveguides with distinct properties and successfully demonstrated 
all the declared capabilities. 

The core of my Thesis begins in Chapter [2] with general discussion of the mathe- 
matical formalism of cross-correlation imaging. Application of the method to charac- 
terization of several optical waveguides is demonstrated in Chapter [3) The summary 
of my experimental and theoretical work is presented in Chapter [4j 
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Chapter 2 

Mathematical formalism of C 2 -imaging 



In this Chapter, I discuss the mathematical formalism of C -imaging (Barankov and 



Ramachandran, 2012), which forms the basis of the experimental method applied 
in Chapter [3] for characterization of several optical waveguides. In Section |2.1[ I 
introduce the general framework of waveguide interferometry necessary for basic ap- 



plications of the method. Then, in Section 2.2, I present a modification of the method 



suitable for studies of polarization properties of waveguide modes. Finally, in Sec- 



tion 2.3, I describe phase-sensitive modification of the method. 



2.1 General description of C 2 -imaging 

In this section, I provide the basic description of the mathematical formalism of 
C 2 -imaging method. Some elements of the formalism were already presented in 



Ref. (Schimpf et al., 2011); here, I expand and further develop this early discussion. 

In multimoded optical waveguides, a beam of light propagates as a superposition of 
discrete modes characterized by different propagation constants and various intensity 
patterns of the modes, depending on the boundary conditions. Direct imaging of the 
beam provides the intensity distribution and also the output power averaged over 
superposition of all the modes. The goal of this section is to demonstrate a novel 
interferometric approach that allows characterization of modal powers and intensity 
distribution of every waveguide mode contributing to the superposition state. 

The basic idea of the method is to study the interference of the beam radiated from 



an optical waveguide with an external reference beam and detect different waveguide 
modes in the time-domain by changing the relative optical paths of the two beams. 

2.1.1 Interferometry of optical beams 

In this work, we employ the standard Mach-Zehnder interferometer configuration 
shown in Figure |2-1| to study the interference of the reference and the test beam 
radiated from the signal waveguide. The beam of light from an LED source is divided 
into two beams at the beam-splitter and then, after propagating in the reference and 
test waveguides, recombined at the beam combiner. The beam of light, radiated from 
a single-moded reference fiber, is collimated, while the signal beam of the test fiber 
is focused at the camera. The camera records a stack of images at different positions 
of the delay stage. 




Figure 21: C 2 -imaging setup: LED - light-emitting diode, BS - 
beam-splitter, BC - beam combiner, Delay - delay stage. 
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The electric field at the image plane of the camera is a superposition of two fields 



E(r,t) = E r (r,t) + E s (r,t + r), (2.1) 

where E r (r,£) and E s (r, t) are the electric fields of the reference and signal beams, 
correspondingly. 

The signal beam is delayed by time 

r = d/c (2.2) 

with respect to the reference beam, where d is the free-space difference between the 
two optical paths, and c is the speed of light in vacuum. The time delay is introduced 



by the delay stage shown in Fig. |21| 

The camera records the intensity of light averaged over the exposure time, which 
allows to employ the two equivalent representations of the electric fields in time and 
frequency domains, E(r,cj) = J d£e zc ^E(r,t), according to the well-known identity 
(Parseval's theorem): 

+AT/2 +oo 

I(r,r)= J d,|E(r,,)| 2 =|^|E(r,a;)| 2 , (2.3) 

-AT/2 -oo 

where r = (x, y) defines the position of the camera pixel in the imaging plane, and AT 
is the exposure time of the camera, assumed to be much larger than the characteristic 
period of the electric field oscillations. 



Upon substitution of the electric field (2.1) into Eq. (2.3), we arrive at the expres 



sion for the average intensity that contains the background intensities of the reference 
and the signal beams Iq : independent of the time-delay r, and also the interferometric 
intensity I int 

/(r,r)=/ (r) + / int (r,r), (2.4) 
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with two terms written explicitly as 

+00 

/o(r)=/^(|E r (r,o;)| 2 + |E s (r,c)| 2 ), (2.5) 

— CO 

OO 

I int (r,r) = 2Re J dt E* r (r, t)E s (r, t + r) 



-OO 



2Re / ^E;(r,o;)E s (r,a;)e-^, (2.6) 



According to these equations, the interferometric signal Ii nt ~ E*E S ~ e lLf>rs IUJT is 
sensitive to the optical path difference cp rs — (p r — tp a and also the time-delay r. 



2.1.2 Analysis of interferometric signal 



The electric fields of the reference and signal beams in Eqs. (2.5) and (2.6) are given 
by 

E r (r, u) = e r (r,u)A r (uj)e i ' Pr , 

E s (r,u) = J2 a m e m(r,ou)A m (u)e ilpm , (2.7) 

m 

where the electric field E s is a superposition of waveguide modes E m with real modal 
amplitudes a m . The power of the m-th mode is given by 

p m = a 2 m . (2.8) 

The spectral properties of the reference beam and waveguide modes are encoded in 
the corresponding functions A r (cj) and A m (lj). 

Every mode of the signal beam and the reference beam is characterized by the 
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intensity distribution: 



I re f(r,u) = |e r (r,cj)| 2 , J m (r,o;) = |e m (r,cj)| 2 . (2.9) 

The intensity distribution 7 re /(r, uj) of the reference beam is imaged directly by block- 
ing the signal path of the interferometer. A similar direct measurement of the signal 
beam, done by blocking the reference path, produces only the total intensity distri- 
bution I tot = I m of all the modes propagating in the signal beam, and, thus, does 
not provide access to the intensities of individual modes. 

The goal of this section is to demonstrate that C 2 -imaging allows to separately 
measure the intensity distribution and relative power level of every mode of the signal 
beam, employing the property of the modes to experience different phase shifts upon 
propagation in a waveguide. 

Specifically, the m-th mode, after propagation in the test fiber of length L, acquires 
the phase 

<p m = l3 m (u>)L, (2.10) 

where f3 m (uj) is the propagation constant of the mode. 

The phase of the reference field, propagating in a single-mode reference fiber of 
length L r , is given by 

(p r = f3 r (uj) L r , (2.11) 

where f3 r (uj) is the propagation constant, and L r is the length of the reference fiber. 

The transverse phase-profile of the reference beam is assumed flat at the imaging 
plane, which is satisfied by the collimation of the beam and also ensuring that the 
imaging plane coincides with the neck of the collimated beam. Regarding the signal 
beam, the near-field image of the output facet of the signal fiber is recorded, by 
focusing the beam at the imaging plane. 



By substituting the electric field of Eq. (2.7) into Eq. (2.6), we arrive at the 



expression for the intensity /^(r, T )) which accounts for the interference between the 
reference field and the individual modes: 

+00 

J in i(r,r) = 2Re^ J ^(r.wKH^ir.wj^He^-^* (2.12) 



-OO 



The analysis of this expression simplifies if one makes several reasonable assumptions. 
In particular, since we are mainly interested in optical measurements, typically per- 
formed with relatively narrow spectral widths around the central frequency of the 
source, it is usually safe to assume that the electric fields do not significantly vary in 
this spectral window 

e(r,cj) « e(r,cj ) , (2.13) 
where ujq is the central frequency of the light source. 



Then Eq. ( |2.12D simplifies: 

I in t(r,r) = 2Re ^ a m e^(r, u )e m (r, uo)G rm (r + r rm ) e ~ l0rm , (2.14) 

m 

where I introduced the joint coherence function 

G rm (t) J 2^ c S rm {yt) , 

S rm (Q) = A* r (n)A m (^)e- i(frm{n) , (2.15) 

and defined the spectral functions with respect to a shifted frequency 

n = uj-u . (2.16) 

The joint coherence function depends on the differential group-delay of every mode 
propagating in the signal beam with respect to the reference beam 

T mr L j Vg rm L r jvg rr . (2.17) 
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The group-velocity dispersion mismatch between the two beams is accounted for 
by the frequency-dependent phase 



(p. 



mr 



J2(€ ] L-/3i k) L r )n k /k\, (2.18) 



k>2 



where ft(® stand for the Taylor-coefficients of the mode-propagation constant ft cal- 
culated at the central frequency u)q of the source. 

The dominant contribution to the time-delay variation of the interferometric in- 
tensity comes from the zero-order optical path mismatch of the beams defined at the 
central frequency of the source 

6 mr = /3^L-/3^L r -u T. (2.19) 

The joint coherence function describes the impact of the frequency-dependent 
mode-propagation constant and the spectral properties of the fields on the interfero- 
metric signal. In particular, the group-delay r mr defines the position of m-th mode 
in the interferometric trace recorded by the camera at a given transverse position 
in the imaging plane, and the shape of the trace is determined by the joint spectral 
function of the reference and the signal beams S rm . In cases when waveguide modes 
do not undergo spectral filtering (e.g. due to a long period grating embedded in the 
waveguide), the spectral functions of the two beams are identical, A = A m = A r , 
and the joint spectral function is given by S rm = Soe~ lcprrn , where So = \A\ 2 is the 
spectrum of the light source. 



The general expressions (2.14) and (2.15), describing interference of the signal and 



the reference beams, form the basis of further theoretical analysis and experiments. 
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2.1.3 Determination of mode dispersion and relative group delays 

The structure of the coherence function defines dispersive properties of the modes. 
Specifically, the complex- valued joint coherence function may be represented as a 
product of the real- valued amplitude and the phase factor: 

G mr (r) = Onr(r)e-^M (2 . 2 0) 

which certainly depend on the detailed structure of the source spectrum. In general, 
however, these are slowly varying functions on the time scale of the electric field 
oscillations 2n /u)q. 

An important property of the joint coherence function is the normalization, which 
is independent of the dispersive properties of the modes. Indeed, it is determined 
only by the spectral properties, as it follows from the identities 

oo oo oo 

I dr\G mr (r)\ 2 = J drGl r {r)= f ^ \A m {Sl)\ 2 |A r (f2)| 2 . (2.21) 

— oo — oo — oo 

Normalization condition also indicates that the amplitude of the joint coherence 
function Gmr(t) has finite extent A m in the time domain defined as 

t £ A m : \t\ < 5t m ~ 1/Acj, (2.22) 

which is dictated by the smallest spectral width of the reference and the signal beams 
Auj. The time-extent 5t m is usually mode-dependent, due to the dispersion effects. 
The separation of time-scales, 5t m 3> 1/ujq suggests time-averaging of the square of 



the interference term (2.14) on the time-scale 1/ljq <C At <C St m . After the averaging, 
we obtain 

r+At/2 

Pr(r,r) = ^ J dtJ?a(r,t)/J r e/(0 = ^Pm|^e m | 2 ^ r (r - r mr )/ m (r), (2.23) 

r-At/2 m 
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where p m is the modal power, e r and e m are the polarization unit vectors of the 
reference and signal beams, and 

7 re /(r) = |e r (r,cj )| 2 , / m (r) = |e m (r,cj )| 2 (2.24) 

are the intensity distributions of the reference and signal beam at the central frequency 
of the source. 

By measuring the interference signal for a full set of polarization states of the 
reference beam, one finds the intensity as a function of the time-delay: 

V(v,t) = £> r (r,r) = Y,P™GL(t ~ r mr )I m (r)- (2-25) 

r m 

This expression indicates a strategy for separating waveguide modes in the time- 
domain. In particular, provided the differential group delay of the two nearby modes 
is larger than the characteristic extent of the coherence function, i.e. r m ^ m+ i > 
the interference peaks corresponding to these modes are well-separated in the time- 
domain. In this case, the intensity P(r,r) at every position r in the imaging plane 
exhibits a series of well-separated peaks located at r = r mr , having the shapes deter- 
mined by the spectral functions of the modes. 

Therefore, by employing an appropriate fitting procedure of the joint coherence 
function to the measured interferometric signal, one finds the relative group delays 
and determines dispersive properties of the modes. In cases when the dispersion of 
the reference beam is unknown, one can measure the interferometric signal at several 
different lengths of the reference waveguide, and then use the resulting set of equations 
to determine dispersion characteristics of all the modes, including the reference mode. 
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2.1.4 Reconstruction of modal weights and intensity distributions 

The structure of the joint coherence function allows determination of the modal 
weights without detailed knowledge of the source spectrum. Indeed, the modal inten- 



sity distribution and power of every mode can be found by first, integrating Eq. (2.25) 
over the extent A m and then integrating over the imaging plane: 

/m(r) = J rfrP(r,r), 

Am 

Pm = I dr [ dr7>(r,r), (2.26) 



where it is assumed that the modal intensity and spectral function (same for all 
modes) are normalized to unity, i.e. 

J ^\S mr (ty\ 2 = 1, J dvl m {v) = 1. (2.27) 

The relative powers p m of all waveguide modes can be reliably extracted from the 
interference signal only when the coherence peaks do not significantly overlap. When 
some peaks do overlap, e.g. due to their relatively small differential group delay, the 
outlined procedure may still be used. However, in this case it will provide the relative 
power of the corresponding group of overlapping peaks. The formalism presented 
above can not discriminate the degenerate modes. Below, I introduce a polarization- 
sensitive modification of the method that allows the reconstruction of degenerate 
modes, using their polarization properties. One example of such reconstruction is 
demonstrated in Chapter |3j Other modifications, accounting for mode-dependent 
spectral properties, are also possible. 
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2.1.5 Gaussian model 

The general description of the beam interferometry presented in the previous section 
is illustrated in this part, assuming same spectral characteristics of the signal and 
reference beams, A = A m = A r , and using a gaussian model of the light source 

S(Q) = S exp [-(n/AQ) 2 ] , Q = uj - uj , (2.28) 

centered at ujo with the spectral width All To simplify the calculation, I consider 
only the effects of group- velocity dispersion 

¥W« (/3^L-/3^L r )n 2 /2. (2.29) 



Under these conditions, the integral in Eq.(2.15) can be analytically calculated, 



which leads to the following expression for the interferometric intensity 



V(r,r) = ^p m J m (r)— = exp 



(r-r mr ) 2 AQ 2 



(2.30) 



2(1 + ^L) 

where the dispersion effects are included into the dimensionless parameter 

d mr = (f3$L - /3^L r )AQ 2 /2, (2.31) 
and the modal intensities I m {v) are normalized. 



Eq. (2.30) indicates that relative height of the interference peaks strongly depends 
on the dispersion effects: for larger dispersion, the time extent of the peak increases 
while its height correspondingly decreases, so that the normalization of the coherence 
function is preserved. This observation is generally true for any spectrum properties 
of the reference and signal beams and suggest the use of the integral characteristics 
instead of the peak value for measurements of the modal powers. 
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2.1.6 Dispersion compensation 

The temporal extent of the interference peaks corresponding to different modes defines 
the temporal resolution of C 2 -imaging. Most clearly this limitation appears in the 
analysis of the gaussian model, where the temporal full-width at half-maximum of 
the interference peak 

&t fwhm = V8h^2^/l + d 2 mr /An. (2.32) 

is determined by the spectral width Aft of the source and group-velocity dispersion 
mismatch between the two beams d mr . Obviously, this expression provides an esti- 
mate for the temporal resolution of the method, assuming that a pair of the modes 
with small relative group delay have similar group- velocity dispersions. 



As it follows from Eq. (2.32), temporal resolution approaches the limiting value 

At fwhm = V8 In 2/ Aft, (2.33) 

when the lengths of the reference and the test fibers are matched, i.e. 

L opt = L^I^\ (2.34) 

so that the group- velocity dispersion is fully compensated. 

For an arbitrary non-gaussian spectrum, the condition of the dispersion-compensation 
remains the same. Indeed, the dispersion effects are accounted for by the frequency- 
dependent phase 

¥W = E W L - Pr k) L r ) n k /kl, (2.35) 

k>2 

The structure of the phase shows that the impact of dispersion is minimized when 
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the length of the reference waveguide satisfies 



J opt 



Lffilfi k \ (2-36) 



where k is the order of the first non-vanishing dispersion term. The dispersion effects 
usually become noticeable already at k — 2, and, as expected, this condition then 



reproduces Eq. (2.34). It's worth noting that conditions (2.34), and (2.36) can be 
satisfied only when the ratios of the corresponding dispersive coefficients fiffl //3r > 0. 

In practice, it is possible to compensate the dispersion effects for a single mode or 
a group of modes with similar dispersion characteristics. The residual dispersion of 
other non-compensated modes results in their temporal broadening. 



Certainly, the limit of the temporal resolution (2.33) can only be achieved for 
purely parabolic dispersion. In more general situations, higher-order dispersive terms 
may become important. For example, when the group- velocity dispersion is fully 
compensated, the third-order phase term contributes a correction to the spectrally 
limited temporal resolution of the order 

dt\ * LAQ 3 - /W) / (6/5( 2 )) , (2.37) 

which depends on the spectral properties of the system and the length of the test 
waveguide. Under usual conditions, this correction is expected to be small. 

2.2 Polarization-sensitive imaging 

In this section, I describe how C 2 -imaging can be used to reconstruct polarization 



state of the modes propagating in few-moded optical waveguides (Barankov et al., 



2012b). Specifically, I show that the measurement of the cross-correlation signal 
for six different polarization states of the reference beam allows full characterization 
of the polarization state of every mode of the test beam in terms of the spatially- 
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dependent Stokes parameters. A specific implementation and application of this 
method is demonstrated in Chapter [3j 

2.2.1 Polarization states 

In spinor notations, every polarization state of m-th mode can be expressed as a 
vector with two components 

=**{\E&)- (2 ' 38 » 

where (p — a y — a x is the relative phase of the complex amplitudes E x = \E x \e lOLx and 
E y = \E y \e tay . The absolute phase, without loss of generality, can be neglected. 

The complete characterization of the polarization state requires identification of 
the following three pairs of orthonormal basis states. The first pair is formed by the 
states linearly polarized along x- and y-axis: 

I#>=(J). I V )=(?)- ( 2 - 39 ) 
Another pair is obtained by rotating the basis by 7r/4 radians: 

| + > = 4m. I") = 4f \\ (2.40) 



V2 v 1 J ' V2 v - 1 

The last pair is given by the left and right circularly polarized states: 

i£> ^(-0- w= ^(0- (2 - 4i) 

Experimentally, projection onto these six states can be realized using combinations 
of a linear polarizer, half-wave and also quarter-wave plate, characterized by the 
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following matrices 



T LP {6) 
Thwp(9) 
Tqwp(0) 



cos 2 9 sin 9 cos 9 
sin 9 cos 9 sin 2 9 

- cos 20 - sin 29 

- sin 29 cos 20 

z cos 2 9 + sin 2 (z — 1) sin 9 cos 9 
(i — 1) sin 9 cos z sin 2 9 + cos 2 



(2.42) 



where 9 determines the angle of the fast axis of the corresponding polarizing element. 

A specific implementation of the projections onto these states is demonstrated in 
Chapter [3j 

2.2.2 Stokes parameters 

An arbitrary state of polarized or unpolarized light is characterized by the Stokes 
vector S = (So, Si, S2, S3) with the following components: 



So 
-Si 
S 2 
S3 



(l^| 2 ) + (|^| 2 ), 
(l^| 2 )-(|^| 2 ), 

2(\E x E y \cOS if), 

2(\E x E y \ sin ip), 



(2.43) 



where averaging (...) is taken with respect to an ensemble of measurements. When 
E x , E y and ip do not vary within the ensemble, the light is completely polarized. 

2.2.3 Polarization distribution of waveguide modes 

In C 2 -imaging, the intensity distribution function is proportional to the scalar product 
of the polarizations of the reference beam e r and the m-th mode e m (r) = (E x (r), E y (r)) rj 
(transposed vector) of the signal beam: 



V rm (r,r) = p m I m (r) \e* r e m (r)\ 2 g^ m (r - r mr ), 



(2.44) 
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The intensity distribution J m (r) is obtained from this expression by summation over 
a complete set of reference states (e.g. \H) and and then integrating over the 
temporal extent of the mode A m , as described earlier in this Chapter. 

The polarization state of every mode, corresponding to a peak value of the coher- 
ence function Q(r — r mr ) at r = r mr , can be reconstructed from the measurement of six 
different polarization states of the reference beam e r . The number of measurements 
is reduced to four when one accounts for the symmetry of the basis states. 

It is convenient to introduce polarization distribution function 

U rm (r) = J rfrP rm (r,r)/(p m / m (r)) = |e;e m (r)| 2 , (2.45) 

Am 

sensitive to the overlap of the reference polarization and spatially dependent mode- 
polarization. 

Straightforward calculations carried out for six spatially uniform polarization 
states of the reference beam lead to the following set of equations: 



U H = 


\E X 2 , 


U v 




\ E y\ 2 ^ 




u + = 


\(\E, 


r+ 


\ E y 


2 ) + \ E X E y\ 


COS (£, 


U- = 




r+ 


\ E y 


2 ) — \E x E y \ 


COS (f 


u L = 


lm 


r+ 


\ E y 


2 ) — \E x E y \ 


sin 99, 


u R = 


\(\E, 


r+ 


\ E y 


2 ) + \E x E y \ 


sin 99, 



with the obvious symmetry identities Uo = Uh + Uv = U+ + U- = Ul + Ur. 

The spatial distribution of the Stokes vector of the ra-th mode is directly obtained 
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from these relations, by averaging over the ensemble of measurements: 



S (r) = {U H ) + {U V ) = {U ), 

Si(r) = {U H )-(U V ) = 2{U H )-(U ), 

S 2 (r) = (U+) - = 2(U+) - (Uo), 

S 3 (r) = {U R )-{U L ) = 2{U R )-(U ) (2.47) 

In the case of a polarized light, the electric fields do not fluctuate within the measure- 
ment ensemble, and the formalism simplifies due to the identity Sq = + + £f 
limiting the polarization vector to the surface of the Poincare sphere. 

This formalism has been used to obtain the spatial distribution of the Stokes 
vector in the case of vortex states, as discussed in Chapter [3) 

2.3 Phase-sensitive imaging 

Spatial phase (p m (r) is one of the key characteristics of waveguide modes, illustrating 
the complex nature of electric field: 



e m (r,cj ) = V im(r) exp [-i<p m (r)] e m (r), (2.48) 

Although in some trivial cases an irrelevant factor, the phase is essential for under- 
standing the vortex modes recently obtained in optical waveguides. 

In this section, I describe how one can use phase-sensitive C 2 -imaging for recon- 



struction of the spatial phase profile of waveguide modes QBarankov et al., 2012a). In 



particular, the intensity distribution of the interferometric signal contains overlap of 
the reference and 771-th waveguide modes: 

Iintfar) = 2Re ^ a m e;(r, o; )e m (r, o; )G m (r - r mr ) e"^, (2.49) 
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The intensity distributions of the reference beam and the m-th mode, as well as 
the polarization distributions can be obtained by the methods described earlier in this 



Chapter. It becomes possible, therefore, as it follows from Eq. (2.49), to introduce 
the phase function that accounts only for the impact of the temporal dependence of 
the coherence function and the spatially dependent phase of the m-mode: 

X rm (r, r) = Grm(r) cos [cp m (r) - uj t + i/; rm (r) + rm (r mr )] , (2.50) 

where, to simplify the notations, time variable was redefined: r — r mr — >■ r. Here, 
^ m (r) and Q rrn {^) are the phase and magnitude of the coherence function G> m (r), 
varying on the time-scale significantly larger than the period 27t/ujq of the electric 
field. 



As it follows from Eq. (2.50), one can use the following strategy to extract the 



spatial phase (p m (r) from the phase function (2.50). For a given position r#, we 



identify moments of time r max and r , such that r max — r < tt/(2ujo) and 

-^rm(^ '* i T~max) ~ QrvrtK^vnax) ) -^771(^**5^0) = . T max Tq ^ 7I~/(2cJo)- (2.51) 

As a result, one finds two complimentary functions 

2™(r, r max ) = g r m(r max ) cos [(p(r) - <p(r*)] , 

2rm(r, r ) = Grm(ro) sin [cp(r) - <p(r*)] , (2.52) 

which allow reconstruction of the spatial phase 

(p(r) - (p(r*) = -i In [X rm (r, T max ) /Q rm {T max ) + iX rm (r, r )/^ rm (r )] . (2.53) 

This approach has been used for reconstruction of spatial phases of the vortex 
modes, as demonstrated in Chapter [3| 
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Chapter 3 

C 2 -imaging in experiments 

In this Chapter, I discuss application of C 2 -imaging for characterization of several op- 



tical waveguides (Barankov and Ramachandran, 2012). In Section 3.1, I demonstrate 



the consistency of the method in retrieving modal content of a few-moded specialty 
fiber with an embedded long-period Bragg grating, acting as a mode-converter. In 
this system, one waveguide mode is converted into another as a function of wave- 
length, by the resonant mode-coupling mechanism. Thus, an independent verifica- 
tion of the modal content measured using C 2 -imaging method becomes possible by 
an independent spectroscopic measurement of the mode- conversion efficiency of the 
mode-converter. 



Next, in Section 3.2, the dispersion-compensation is discussed and applied for 



analysis of the modal content of large-mode area fibers. 



Polarization-sensitive modification of the method is applied in Section 3J3 for 
reconstruction of mode-specific polarization distribution of vector modes propagating 
in a specialty fiber. Specifically, the Stokes parameters of the modes are compared to 
the theoretical predictions, providing the first to-date verification of the polarization 
patterns of these modes in multimoded fibers. 

The specialty fiber also supports two degenerate vector modes carrying orbital an- 
gular momenta resulting in the characteristic discontinuity of their spatial phase dis- 



tribution. In Section 3.4, I apply phase-sensitive modification of C -imaging method 



for complete characterization of the vortex modes, and obtain their relative weights, 



polarization and phase properties. 

In Section 3.5, I employ C 2 -imaging for identifying the mechanism of resonant 



mode-coupling observed in coiled large-mode- area leakage- channel fibers. 

Regarding the experimental setup, the modal content of higher-order-mode fiber 
with embedded long-period grating, and also the measurements of large-mode-area 
fibers illustrating the concept of dispersion-compensation were carried out using the 
basic experimental setup discussed in Chapter [2j The polarization and phase-sensitive 
measurements were carried out using a later-developed polarization-sensitive modifi- 



cation of the setup discussed in Sections |3.3| and |3.4 



The content of Sections |3.1|and 3.2, with some clarifying textual changes, closely 



reproduces the corresponding parts of the journal article (Schimpf et al., 2011), to 
which I have been a key contributor. D. Schimpf and I have contributed equally to 
work reported in these sections. 

3.1 Modal content of higher-order mode fiber with long- 
period grating 

The specialty higher-order mode (HOM) fiber supplemented by a turn-around-point 
long-period grating (TAP-LPG) unit allows direct comparison of modal content re- 
trieved by C 2 -imaging method to independent spectroscopic measurements of the 
conversion efficiency. Since the higher-order modes supported by this fiber demon- 
strate distinct dispersive behavior, the same fiber is used to illustrate the impact of 
dispersion on the interferometric signal. 

In the first set of experiments we employed free-space path for the reference beam 
(the dispersion of optical elements positioned in this path was negligible) . A schematic 
of the setup is shown in Fig. |3T| The few-mode fiber under test is the final element of 
a module (L=0.6 m) consisting of a single-mode fiber, a TAP-LPG, and the higher- 



23 



delay dispersion beam- band- 
stage compensating combiner pass camera 
<<— ► SM fiber 



mode 1 



mode 2 



SLD 




Figure 31: Schematic of the experimental setup (SLD: superlumines- 
cent diode), and illustration of the cross-correlation trace expected at 
one pixel of the stack of images. 



order mode fiber (L=0.4 m) flRamachandran, 2005 ; Jespersen et al., 2010). In this 
fiber, the LPG spectrum characterizes the mode conversion efficiency from the LPqi 
core-mode to the LP 2 core-mode, which is used as an independent reference for 
comparison against the ratio of modal weights obtained by C 2 -imaging method. 

The source spectrum, provided by an LED, was filtered by a bandpass with the 
central wavelength of about 780 nm and bandwidth of about 4 nm. Since the LPG 
mode-conversion bandwidth exceeds 20 nm, the wavelength-dependent LPG spectrum 
has a negligible impact on the spectra of HOM's. 

Figure |3-2| (a) shows an example of cross-correlation trace between the reference 
field and the output of the fiber, integrated over all pixels of the camera. The peaks 
in the trace correspond to the two different modes in the HOM fiber. We analyze the 
data using elements of the theoretical analysis discussed in Chapter [2j 

The envelope of the cross-correlation trace is obtained from the stack of images 



recorded by the camera, as shown in red color in Fig. 3-2 (a). Then, the amplitude of 
the coherence function Q vm {r) is fitted to the extracted envelope, separately for the 
two peaks, as Figs. |3-2[b) and (c) demonstrate. 



The difference in shape of the two interference peaks reflects distinct modal disper- 
sion of the corresponding modes. The side-lobes around the dominant peaks (clearly 
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Figure 3-2: (a): Cross-correlation trace for the entire image (data is 
offset corrected) for the bandpass centered at A center =780 nm. (b) and 
(c): fit of the model to the envelope of the experimental data, for the 
first and second peak, corresponding to LP i and LP 2, respectively. 

visible for the LP 02 mode) are due to the steep spectral edges of the filtered spec- 
trum. By fitting Q rrn {j) to these data, we obtain the relative group-delays and the 
dispersion of the two modes. 

The dependence of the group-delay and dispersion on the wavelength is measured 
by shifting the central wavelength of the filtered spectrum (achieved by changing the 
incidence angle of light on the bandpass). The results are shown in Figs. |3-3| (a) and 
(b). For comparison, we have also simulated the modal properties of the tested fiber 
using a scalar mode-solver. 

Figures |3-4[ a) and (b) show the reconstructed intensity distributions of LP i and 
LP 2 modes that match the expected pattern of the modes. 

The relative (dispersion-corrected) weights of the normalized modes are charac- 
terized using the multi-path interference (MPI) value 

MPI = 101og 10 (p m /p L p 01 ) , (3.1) 
a useful measure of the relative power of a higher-order mode, p m to that of the 
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Figure 3-3: (a) Group-delays, and (b) Dispersion values of the two 
modes as a function of center wavelength of the bandpass. 
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Figure 3-4: (a) and (b), reconstructed LP i and LP 2-mode (gamma- 
adjusted) at a center wavelength of 780 nm, (c) multi-path interference 
(MPI) values as a function of center wavelength of the bandpass. 



fundamental mode p rp . 

We demonstrate the accuracy of the C 2 imaging in Fig. |3-4| (c), where the multi- 
path interference (MPI) value is plotted as a function of the wavelength. The obtained 
MPI values match very well with the mode conversion efficiencies independently mea- 
sured by recording the LPG spectrum. 
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3.2 Dispersion compensation for imaging of large-mode area 
fibers 



The advantages of dispersion compensation are illustrated in characterization of large- 
mode-area (LMA) fibers supporting modes with similar dispersive characteristics. As 
I argued in Chapter [2j in this case, it is possible to achieve high temporal resolution 
of small intermodal group delays, although, obviously, at the expense of spectral 
resolution. 

The temporal resolution of C 2 -imaging is related to the width of the coherence 
function G mr (r), which is determined by the spectral properties of the light source. 
For a given shape of the spectrum, the temporal resolution Atfwhm (defined here 
as the FWHM of Q mr {r)) can be calculated as a function of the bandwidth of the 
spectrum and the residual group- velocity dispersion (GVD) Aip( 2 \ 
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Figure 3-5: (a) Temporal resolution as a function of the FWHM 
spectral bandwidth of a Gaussian spectrum (for GVD value of tp^ = 
O.lps 2 ), (b) and as a function of both FWHM bandwidth and GVD. 



We modeled the dependence of the temporal resolution for a Gaussian spectrum 
centered at 1060 nm and residual GVD of Atp^ = O.lps 2 (corresponding to about 
4 m of LMA fiber and in the absence of the dispersion compensating fiber in the 
reference path). The results are shown in Fig. |3-5[a). 
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For small band widths, dispersion plays a minor role, and the resolution is gov- 
erned by the coherence time, which can be defined as Ar^°^ M = 8ln(2) / Aljfwhmi 
where Aujfwhm is the width of the spectrum at the angular frequency (Aujfwhm ~ 
27tc AX fwhm /\1). 

In the absence of dispersion, broad spectra may be used to achieve maximal tem- 
poral resolution. However, when the interferometer is not dispersion balanced, the 
broadening of cross-correlation signal is dominated by the dispersion, i.e. At^ hm = 
fwhm • Nonetheless, even in this case, it is possible to identify an optimal 
bandwidth of the spectrum (i.e. appropriate choice of the bandpass) corresponding 



to the optimal resolution, as Fig. 3-5 (b) clearly illustrates. The horizontal line in Fig. 
3-5|(b) highlights the parameter configuration shown in Fig. |3-5|(a). 



The temporal resolution may be noticeably improved in a dispersion-matched 



interferometer, a well-known result in optical coherence tomography (OCT) (Wo- 



jtkowski et al., 2004). In our setup, we realize this idea by employing a single-mode 
fiber in the reference path, which balances the dispersion of the test fiber. Certainly, 
the dispersion can be exactly matched for a single mode or group of modes with 
similar dispersive characteristics. For the other modes, the residual dispersive phases 
cause broadening of the corresponding coherence peaks in the cross-correlation trace. 

The shape and smoothness of the spectrum also have an impact on the cross- 
correlation trace, and thus, affect the temporal resolution. These effects have also 



been discussed in the context of OCT (de Boer et al., 2001). 

For characterization of large-mode area (LMA) fibers, in which all the modes have 
similar magnitudes of chromatic dispersion and the relative delays are on a picosecond 
(or less) timescales, the results of dispersion-matching are particularly illuminating, 
as we demonstrate below. 

We demonstrate dispersion compensation for characterization of the polarization- 
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maintaining LMA fiber with a core diameter of approximately 27.5 \±m and a numeri- 
cal aperture (NA) of 0.062, having the length of 5 m. We use a polarizer and a pair of 
half- wave plates to ensure launch of the beam into one of the dominant polarization 
states of the fiber. The dispersion is matched by using 4.08 meters of single-mode HI- 
1060 fiber in the reference path. This length has been determined by cutting the input 
of the reference fiber until the width of the dominant peak in the cross-correlation 
trace matched the coherence function calculated for the full spectrum of the source. 
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Figure 3-6: (a): Full spectrum of the source and the 5-nm bandpass, 
(b): the envelopes of the cross-correlation traces. 
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Figure 3-7: Reconstructed mode profiles in the order of temporal 



delays shown in the cross-correlation trace of Fig. 3-6 (b) for the case 
of the full spectrum of the source. 



The effect of dispersion-matching is clearly visible when comparing the cross- 
correlation traces, obtained using the full spectrum of the source in Fig. |3-6|(a), and 



those recorded using a spectrum filtered by a 5-nm bandpass filter in Fig. 3-6 (b). Most 
importantly, dispersion compensation combined with the use of full spectrum reveals 
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Figure 3-8: (a-c) Output near-field images of the tested fiber, obtained 
for different in-coupling conditions, and corresponding changes in the 
cross-correlation trace. 

the pairs of LPu and LP 2 i modes, which overlap in the traces obtained using the 
filtered spectrum. The corresponding temporal resolution approaches values smaller 
than 300 fs. The impact of the shape of the spectrum on the cross-correlation trace 
is also noticeable: since the filtered spectrum has steep edges, the corresponding 
cross-correlation trace shows characteristic ringing (especially around the peak of the 
LP i- m °de). In contrast, the smooth full spectrum results in a cross-correlation trace 
without spectral artifacts. 

Fig. |3-7| demonstrates intensity distributions of all modes, arranged in the order 
of their time delay in Fig. 3-6| (b). The dispersion-corrected MPI values of the modes 
shown in Fig. |3-6| (b) (full spectrum) are estimated as -14.6 dB and -16.6 dB for the 
two LPn modes, and -29 dB and -19 dB for the two LP 2 i modes. Interestingly, 
the dimensionless V-parameter of the test fiber is about 5, so the next higher LP 2- 
mode should also be present in the trace, but it is not observed in our measurements. 
Coiling of the fiber to a diameter of around 30 cm may be responsible for stripping 
off of the mode. 

It is instructive to study the change of the modal pattern for different in-coupling 
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conditions. A qualitative observation based upon the image of the output of the fiber 
in Fig. 3-8[ a), obtained when the reference beam is blocked, suggests fundamental- 



mode operation. Conversely, C 2 -imaging reveals that in fact the beam contains a 
significant amount of power in the higher-order modes. The evolution of the modal 
content for other in-coupling conditions measured by C 2 -imaging is illustrated in 



Fig. 3-8, Here, the near-field image of the fiber-output in Fig. |3-8[ b) indicates a 
stronger excitation of the "slow" LPn-mode, while Fig. |3-8[ c) shows a stronger "fast" 
LPn-mode. Interestingly, in these two situations, LP i-mode still carries most of the 
modal power. 

These measurements also reveal that temporal splitting between the LP21 modes is 
more pronounced compared to the splitting between the LP\\ modes. It is an indica- 
tion that modes with higher orbital angular momentum (i.e. LPi m modes with higher 
I) are more susceptible to the birefringence effects in this polarization maintaining 



fiber (Golowich and Ramachandran, 2005). 



3.3 Polarization reconstruction of vector modes 

One of the major advantages of C 2 -imaging method compared to other waveguide 



characterization techniques (iNicholson et al., 2009) is the flexibility in choosing the 



polarization states of reference and signal beams independently of one another, as 
described in Chapter [2} 

In this section, I employ this advantage of the method to characterize polarization 



properties of vector modes TM i, HE 2 \ ) and TE i flBarankov et al., 2012b) prop 



agating in the specialty fiber, where the degeneracy between the modes is lifted by 
the ingenious design of the index profile. The analysis is based on the formalism 
developed in Chapter [2| 

The four vector modes propagating in the fiber can be represented in the following 
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way: 



eTMoi 


= f(r) (x cos ip + y sin (p) , 




e TE 01 


= /(r) (f sin 99 — y sin 99) , 




e+ 


= f(r)e^(x + iy)/V2, 




e_ 


= f(r)e-^(x-iy)/V2, 


(3.2) 



where f(r) is the radial distribution of the modes that has characteristic "donut" 
shape; x and y are the unit vectors in the image plane. It is convenient to employ 
the cylindrical coordinate system, x = r cos cp : y = rsin(^, where r is measured from 
the axis of the fiber, and p is the azimuthal angle. 

The last two modes e± ~ e ±Lp are the degenerate orbital-angular-momentum 
(OAM) states characterized by the spatial phase p that encodes the orbital mo- 
mentum L z — ±1. Due to the spectral degeneracy, an arbitrary linear superposition 
of the modes defines HE 2 \ mode: 

^he 21 = c+e + + c_e_, (3.3) 

where c + and c_ are the complex- valued amplitudes of the modes normalized to unity, 
|c + | 2 + |c_| 2 = 1. The amplitudes depend on the in-coupling conditions to the fiber. 

The polarization-sensitive C 2 -imaging method is implemented using the experi- 
mental setup shown in Fig. |3-9[ I employ a pair of linear polarizers and half-wave 
plates to control the polarization state of the reference beam. The in-coupling con- 
ditions to the test fiber is also determined by a similar set of polarizing elements. 
Polarization selection at the output of the system is provided by a combination of a 
quarter- wave plate and a linear polarizer. 

I used an LED source centered at about A = 1524 nm with spectral width of ^ 
30 nm. The output beam of the test fiber focused at the camera interferes with the col- 
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Figure 3-9: C 2 -imaging setup: LED light-emitting diode, LP-linear 
polarizer, HWP half-wave plate, QWP quarter-wave plate, BS beam- 
splitter, BC beam combiner. 

limated reference beam having well-defined linear polarization. The in-coupling condi- 
tions remain the same throughout our measurements. The length of the polarization- 
maintaining reference fiber is chosen to match the optical paths of the two beams and 
also to compensate for the modal dispersion of the test fiber. The effect of disper- 
sion has been also minimized by using a spectral filter with a central wavelength of 
A « 1549 nm and width Xfwhm ~ 6 ran. The camera detects a cross-correlation sig- 
nal as a function of the group delay monitored by a computer-controlled delay stage. 
In the measurements, the cross-correlation trace has been recorded for three pairs of 
orthogonal polarization states of the reference beam: linear vertical, linear horizontal, 
linear at +7r/4 radians, linear at — 7r/4 radians, left-circular and right-circular polar- 
izations. The envelopes of the cross-correlation traces have been analyzed to study 
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the polarization properties and the relative power of the modes. 

Among the six polarization states used in the experiment, the left-circular and 
right-circular states are particularly important for identifying the topological charge 
of the OAM state. Specifically, as Fig. |3-10| illustrates, TEqi and TM i modes are in- 
sensitive to changes in orientation of the circular polarizer, as dictated by their spatial 
and polarization properties. In contrast, the two OAM states with the orbital angular 
momentum L z = ±1, have the right and left circular polarizations, correspondingly. 
Therefore, one can detect each of these states using a circular polarizer, as shown in 



Fig. |3-10| The peak values of the thus-recorded cross-correlation traces contain the 
relative weights of the two OAMs, and also the weights of the other two modes. 
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Figure 3-10: Relative modal power as a function of group delay: blue 
line right-circular polarization, red line left- circular polarization 



Polarization-sensitive imaging method described in Chapter [2] allows complete 
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TEi 




Figure 311: Reconstruction of Stokes parameters of TE i mode: (a) 
component So, (b) component Si, (c) component S2, (d) component 
S3. Inset: theoretical prediction. 

reconstruction of the Stokes parameters of TM i, TE$\, and HE 2 i modes and direct 
comparison to the theoretical predictions based upon polarization patterns shown in 



Eqs. (3.2) and (3.3). 



Specifically, for TM i the Stokes vector components are 



S (r) = f(r), S 1 (r) = -f(r) cos(2^), S 2 (v) = /(r)sin(2^), S 3 (v) = 0. (3.4) 



Similarly, for TEqi, one obtains 



S (r) = f(r), S 1 (r) = /(r)cos(2^), S 2 (v) = -f(r) sin(2^), S 3 (v) = 0. (3.5) 



In the most interesting case of HE 2 i mode, the Stokes parameters Si and S 2 encode 
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Figure 3-12: Reconstruction of Stokes parameters of TM i mode: (a) 
component So, (b) component Si, (c) component S2, (d) component 
S3. Inset: theoretical prediction. 

the interference term c + cl = |c + c_|e m of the two degenerate OAM states 



S {t) = f{r),S 1 (r) = 2\c + c_\f(r)cos{2<p + a), 
S 2 (r) = -2|c + c_|/(r)sin(2^ + a), S 3 (r) = \c+\~ - |c_ 



(3.6) 



The intensity distribution of each mode, given by the component So of the spatially 
dependent Stokes vector, as expected, demonstrates characteristic donut shape, as 
illustrated in Figs. |3-11| |3-12[ and |3-13[ 

In summary, the reconstructed Stokes parameters for all vector modes compare 
very well with the theoretical predictions. In addition, phase patterns Si and S2 in 



Fig. 3-13 reveal the relative phase of the degenerate OAM states e±. 
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Figure 3-13: Reconstruction of Stokes parameters of HE 2 \ mode: (a) 
component So, (b) component Si, (c) component £2, (d) component 
S3. Inset: theoretical prediction 

3.4 Spatial phase reconstruction of vortex modes 



Polarization properties of OAM states shown in Eqs.(3.2) suggest a simple method to 



measure their modal power. Indeed, by projecting the output beam to the left and 
right circular polarized states, one obtains the cross-correlation traces of the states 
independently of one another. Moreover, since the other two vector states, TM i and 
TEqi : are insensitive to the change of the circular polarization, these states can be 



easily identified and separated in the cross-correlation traces, as shown in Fig. |3- 10 
One of the key characteristics of OAM states 



e + = f(r)e^(x + iy)/V2, 
e_ = f(r)e-^(x-iy)/V2, 



(3.7) 



37 



is the phase factor e ±lLp corresponding to the two eigen-values of orbital momentum 
L z — ±1 along the direction of propagation. 

In this section, I apply the phase-sensitive C 2 -imaging to measure the phase distri- 



bution of the vortex states (Barankov et al., 2012a). First, I project the interferomet- 
ric signal onto the right- and left- circular polarization states to select the two vortex 
states independent of one another. The resulting interferometric trace contains the 
information about the spatial phase that can be retrieved using simple phase-stepping 
algorithm developed in Chapter [2j 

The cross-correlation signal is measured with high temporal resolution correspond- 
ing to a fraction of the wavelength in spatial translation of the delay stage. The phase- 
stepping algorithm is capable of retrieving the phase information even for relatively 
noisy conditions. 



The results of phase reconstruction of the OAM states are shown in Figs. 3-14 



and |3-15] for L z = +1 and L z = — 1 states, correspondingly. The spatial pattern 
reveals clockwise and counterclockwise phase wrapping, which uniquely identifies the 
two states. The position of the phase discontinuity in Figs . |3 • 1 4] and |3 • 1 5] allows direct 
measurement of the relative phase of the states. 

In summary, OAM states propagating in optical waveguides demonstrate distinct 
polarization and spatial patterns. Using polarization and phase-sensitive C 2 -imaging 
methods, I have completely characterized these modes, revealing their unique polar- 
ization and phase singularities, the modal power and also the relative phase of the 
states. 
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Figure 3-14: Spatial phase of OAM state with L z = +1 that shows 
characteristic phase wrapping in the counter-clockwise direction. In- 
sets: (a) Intensity distribution of the mode, (b) theoretical distribution 
of spatial phase 




Figure 3-15: Spatial phase of OAM state with L z = — 1 that shows 
characteristic phase wrapping in the clockwise direction. Insets: (a) 
Intensity distribution of the mode, (b) theoretical distribution of spatial 
phase 
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3.5 Resonant mode coupling in leakage channel fibers 

For an updated description of this experiment, I refer the reader to a recent preprint 
"Resonant Bend Loss in Leakage Channel Fibers" by R. A. Barankov, K. Wei, B. 



Samson, and S. Ramachandran posted at arXiv:1205. 5584^1 [physics. optics]. Below 



is the original write-up as submitted to the library of Boston University. 

In this section, I describe modal characterization of large-mode-area leakage chan- 
nel fibers, which demonstrate dramatic power loss at certain coiling radius. Using 
C 2 -imaging, I experimentally attribute this anomaly to a new physical mechanism of 
resonant mode-coupling QBarankov et al., 2012c[ ). 

Resonantly enhanced leakage-channel fibers (LCFs) with large mode areas are de- 
signed to provide high-power propagation of diffraction-limited beams in high-power 



fiber lasers (Dong et al., 2009). The microstructure of these fibers is tailored to en- 
hance the loss of higher-order modes (HOMs) while maintaining tolerable loss of the 
fundamental mode, resulting in single-mode operation with large field diameters. In 



the existing designs (Dong et al., 2009), the boundary between the light-guiding core 
and the cladding consists of low-index circular rods, so that index-continuity of the 
boundary in the azimuthal direction is broken. As a result, all modes propagating in 
the fiber are coupled to the radiative modes of the cladding. However, the presence of 
gaps at the interface naturally promotes significant differential power-loss dominated 
by HOMs. A careful study of the mechanism and different possible fiber-designs 
indicates that high confinement loss for HOMs and low confinement loss for the fun- 
damental mode can be achieved, leading to effectively single-mode operation of the 



structures (Dong et al., 2009). 



The physical mechanism underlying the LCF design suggests sensitivity of light 
propagation to the coiling conditions of the fibers. Indeed, it is well-known that coiled 
fibers experience power loss resulting from bend-induced coupling of the guided modes 
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to the radiative modes of the cladding (Marcuse, 1976; Wong et al., 2005; Dong et al., 



2006; Saitoh et al., 2011). In this case, the loss is significantly higher for HOMs than 



for the fundamental mode, as dictated by larger modal overlap of HOMs and the 
cladding modes. Moreover, when the cladding modes become quasi-guided in coated 



fibers (Renner, 1992), a resonant behavior is observed (Murakami and Tsuchiya, 



1978). 



The bend-loss becomes significant for small enough coiling radii (Love, 1989). The 



critical radius has been estimated in the case of photonic-crystal fibers (Birks et al., 



1997) as R c ~ A 3 /A 2 in the short- wavelength limit, where A is the wavelength of light, 



and A is the characteristic core size (Birks et al., 1997; Nielsen et al., 2004). This 
estimate should also hold for leakage-channel fibers characterized by similar geometry. 
For example, for the LCFs with core size A ^ 50 /im, significant bend loss is expected 
at R c rsj 10 cm in the 1 /xm spectral range. 

In addition to the power-loss resulting from direct radiative coupling of the core 



modes, bending of fibers induces inter-modal coupling in the core (Russell, 2006; Love 



and Durniak, 2007). The coupling increases dramatically when the effective indices 



of the two modes approach one another as a function of the fiber curvature. The 
crossing occurs at some critical radius, which, in general, depends on the properties 
of the coupled modes. Interestingly, the crossing between the fundamental and the 
lowest HOM is also expected at the critical radius R c corresponding to onset of 



significant bend-loss QRussell, 2006 ). 

In this work, we explore the effect of resonant mode-coupling on the bend-loss 
in a large-mode area LCF and observe an enhanced power loss at a certain coiling 
diameter. The results of power measurements in the 1 /im spectral range are shown 



in Fig. |3-16[ We observe a dramatic decrease of the output light power at a specific 
coiling diameter, while the power recovers to the levels offset by the usual mecha- 
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Figure 3-16: Output power as a function of coiling diameter. Insets 
(a-c): Output mode profiles at different coiling diameters; (d) Cross- 
section of the LCF 

nism of power loss outside the resonance. At small coiling radii, as expected, one 



obtains significant loss of power in accord with the usual bend-loss mechanism (Mar- 
cuse, 1976; Birks et al., 1997). We characterize the observed anomaly by C 2 -imaging 



method (Schimpf et al., 2011). Specifically, we find that for non-critical coiling radii, 



light propagation in the tested LCF is dominated by the fundamental mode with HOM 
extinction below — 25 dB. Thus, at non-critical radii the power-loss can be explained 



by the usual bend-loss mechanisms (Birks et al., 1997; Russell, 2006). In contrast 



at the critical coiling radius, a higher-order mode dominates propagation, indicating 
resonant coupling of the fundamental mode to HOMs, which, by the design of LCFs, 
immediately radiate out of the core. Interestingly, this mechanism is reminiscent of 



resonant mode coupling observed in coated single-mode fibers (Renner, 1992). While 



simple bend-loss measurements shown in Fig. 3-16 reveal this anomalous behavior 
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C 2 -imaging method allows us to quantify this effect. As a result, this precise char- 
acterization method provides critical feedback for future fiber-designs, in which the 
critical radius R c should be defined for specific amplifier packaging constraints. 

The tested LCF of 285 cm-length has a core diameter of 50 fim and cladding di- 
ameter of 400 fim. The two rings of low-index (fluorine-doped) silica regions shown in 
Fig. 3-17[ a) provide the leakage channel. The core, made of silica, is index matched to 
the outer silica glass. A high-index regular acrylate coating applied to the cladding 
ensures stripping out of the cladding modes. The LCF has been designed to have 
negligible HOM content at lengths greater than 3 m. The input end of the fiber was 
spliced to a single-mode fiber to provide same in-coupling conditions throughout the 
experiments. 

The modal content of the LCF was analyzed using C 2 -imaging method developed 
in Chapter [2j The basic idea of the method is to study the interference of the beam 
radiated from an optical waveguide with an external reference beam, and detect 
different waveguide modes in the time-domain by changing the relative optical paths 
of the two beams. 



Figure [3T7| shows a schematic diagram of the experimental setup based on a Mach- 
Zehnder interferometer. We used a superluminescent diode (SLD) source centered at 
A ~ 1050 nm with a spectral width of ~ 30 nm. The output beam of the LCF (focused 
at the imaging plane) interferes with the collimated reference beam radiated from the 
reference fiber. The length of polarization-maintaining reference fiber is chosen to 
compensate for the optical path difference between the two paths and also to reduce 
the effects of group- velocity dispersion of the LCF. The latter is important since 
we employ a light source with a relatively broad spectrum, which, in the absence 
of dispersion compensation, leads to significant dispersion broadening of the cross- 
correlation signal. In particular, the mode-specific resolution of C 2 -imaging method 
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Figure 3-17: C 2 -imaging setup: LED light-emitting diode, LP-linear 
polarizer, HWP half-wave plate, QWP quarter-wave plate, BS beam- 
splitter, BC beam combiner. Inset (a): Cross-section of the LCF 

is defined by the spectral width of the source and also by the dispersion mismatch of 
the reference and the test modes. In large-mode-area LCFs, the material dispersion 
dominates the spectral broadening, which results in similar dispersion properties of 
HOMs. We chose the length of the reference fiber to match the material dispersion 
of LCF and, thus, reduced the dispersive broadening of all the modes simultaneously. 

The cross-correlation signal V{y ) t) as a function of the group delay r and the 
coordinate r in the imaging plane is detected by the camera for different positions 
d = cr (c is the speed of light in vacuum) of the delay stage: 



Here, the summation extends over the modes propagating in LCF, p m is the relative 
modal power of m-th mode (Y^ m Pm = 1), Qmr{^) is the mutual coherence function 
of the reference and test beams, J m (r) is the modal intensity, and r mr is the relative 




(3.8) 



m 
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group delay of the m-th mode with respect to the reference mode. 

Coiling of the fiber affects polarization properties of the test beam. To extract 
the power of every elliptically-polarized mode, we have recorded the cross-correlation 
trace for two orthogonal polarization states of the reference beam. The resulting trace, 



combining the two measurements, is represented in Eq. (3.8). The modal intensity 
distribution I m (r) of every mode was obtained by integrating this expression over the 
time extent of the mode. The relative power of the modes is encoded in the net cross- 



correlation trace obtained by integration of the spatially-dependent trace in Eq. (3.8) 
over the imaging plane position. 

The result of this procedure, applied to the cross-correlation traces recorded at 
the critical coiling radius, is shown in Fig. |3T8 In this figure, the cross-correlation 
peaks identify the modes propagating in the fiber at the corresponding relative group 
delays. The shape of the peaks reflects the corresponding mutual coherence functions, 
while the peak values encode the relative power of the modes. The insets demonstrate 
the intensities of the reconstructed modes. 

The dependence of the output power on the coiling diameter was measured using 



a power meter. The results are shown in Fig. |3T6[ We observe a dramatic decrease 
of the output light power at a specific coiling diameter. The observed critical radius 
R exp « 11 cm is close to the estimate R c ~ 10 cm, suggesting resonant mode-coupling 
as the mechanism responsible for the anomalous power loss. Indeed, at the resonance, 
the output image demonstrates domination of HOMs as shown in Fig. 3T6| (b), while 
the output images in Fig. 3-16| (a) and 3-16[ c), recorded at coiling diameters outside 
the resonance, indicate single-mode operation. 



C -imaging (Schimpf et al., 2011) provides insight into the resonant behavior. The 



envelope of the integrated correlation trace shown in Fig. |3-18| at the critical coiling 
diameter demonstrates two modes (LP i and LPu) propagating with a relative group 
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Figure 3-18: Relative modal power as a function of group delay at 
the resonance (D=22cm). Insets: output image and the images of 
reconstructed modes 

delay of about 0.2 ps/m, with the fundamental mode contributing only about 15% 
of the total power. In contrast, at other coiling diameters, HOMs are suppressed at 



the power level below — 25 dB, as shown in Fig. 3-19 , The strength of the resonance 
depends on the length of the coiled fiber. In a set of similar measurements conducted 
on the LCF of smaller length of about 180 cm, a relatively shallow resonance was 
found at the same coiling diameter. We expect a significantly deeper resonance for 
longer fiber lengths. 

In summary, coiled few-mode fibers experience bend-induced coupling between the 
core modes and power-loss via direct coupling of core modes to the radiative modes 
of the cladding. Using C 2 -imaging, we explore the interplay of these phenomena in 
LCFs. We identify a new resonant power-loss mechanism in these fibers, in which 
higher-order core modes mediate coupling of the fundamental mode to the radiative 
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Figure 3-19: Relative modal power of LP i and LP\\ (peak values) as 
a function of the relative group delay, for different coiling diameters 

modes. The effect becomes evident at a specific coiling diameter, where we observe a 
dramatic decrease of the output power. Outside the resonance, the power recovers to 
the levels determined by the usual bend-loss mechanism. In this work, we quantify this 
anomaly by C 2 -imaging, thus providing critical feedback for the fiber-based amplifier 
designs with certain coiling constraints. 
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Chapter 4 

Conclusions 

In this work, I have developed a novel interferometric method, called C 2 -imaging, 
suitable for complete characterization of waveguide modes. In particular, using this 
method one obtains the modal weights, relative group delays and dispersion properties 
of the modes. Modifications of the method, involving polarization post-selection and 
the analysis of high-frequency oscillations of interferometric signal, reveal polarization 
and phase distributions of all the modes. 

The developed method has been successfully applied to modal characterization of 
waveguides with distinct modal properties. In particular, in a specialty higher-order 
mode fiber, the method allowed accurate reconstruction of the modes having small rel- 
ative group delays and distinct dispersive behavior. The reconstructed modal weights 
agree well with the mode- conversion of the long-period grating spectrum. The multi- 
path interference values down to — 30 dB can be accurately retrieved. The concept 
of dispersion-compensation was demonstrated in the case of large-mode-area fibers, 
where sufficient temporal resolution was achieved to observe temporal birefringent 
splitting between odd and even higher order modes propagating in the polarization 
maintaining fiber. 

Polarization-sensitive modification of C 2 -imaging method was successfully tested 
for modal characterization of a specialty fiber supporting several vector modes. The 
reconstructed polarization patterns, demonstrating non-trivial space dependence, com- 
pare very well with the theoretical predictions. In addition, I have directly observed 



optical orbital angular momentum states propagating in the fiber, by employing cir- 
cular polarization selectivity of the modes. The observed phase patterns of the modes 
are consistent with the theoretical expectations for the corresponding spatially de- 
pendent Stokes parameters. 

Phase-sensitivity of C 2 -imaging was also demonstrated in the experiments with 
the same specialty fiber. I have directly measured, for the first time to the best 
of my knowledge, the phase distribution of the orbital momentum states, using a 
phase-retrieval algorithm based upon analysis of high-frequency oscillations of the 
interferometric cross-correlation signal. As a result, I have completely characterized 
the vortex modes and identified their relative power and the complex relative phase. 

I have also used C 2 -imaging to explore the resonant behavior of power output in 
large-mode-area leakage-channel fibers. In particular, I have identified the resonant 
mode-coupling as the mechanism responsible for the observed dramatic power loss at 
a specific coiling radius. 

The versatility of C 2 imaging, demonstrated in this work, identifies this method as 
a powerful new tool, suitable for waveguide characterization in optical applications. 

I kindly acknowledge support of this project by ARL Grant No. W911NF-06-2- 
0040, ONR Grant Nos. N00014-11-1-0133 and N00014-11-1-0098. 
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